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. . 1 INTRODUCTION 



ABSTRACT 

Recent observations have gathered a considerable sample of high redshift galaxy can- 
didates and determined the evolution of their luminosity function (LF). To interpret 
these findings, we use cosmological SPH simulations including, in addition to standard 
physical processes, a detailed treatment of the Pop III— Pop II transition in early ob- 
jects. The simulated high-z galaxies match remarkably well the amplitude and slope 
of the observed LF in the redshift range 5 < z < 10. The LF shifts towards fainter lu- 
minosities with increasing redshift, while its faint-end slope keeps an almost constant 
value, a » —2. The stellar populations of high- 2 galaxies have ages of 100-300 (40-130) 
Myr at z = 5 (z = 7 — 8), implying an early (z > 9.4) start of their star formation 
activity; the specific star formation rate is almost independent of galactic stellar mass. 
These objects are enriched rapidly with metals and galaxies identified by HST/WFC3 
(Mjjv < —18) show metallicities « 0.1-Z© even at z = 7 — 8. Most of the simulated 
galaxies at z w 7 (noticeably the smallest ones) are virtually dust-free, and none of 
them has an extinction larger than E(B — V) = 0.01. The bulk (50%) of the ionizing 
photons is produced by objects populating the faint-end of the LF (Mjjv < — 16), 
which JWST will resolve up to z — 7.3. PopIII stars continue to form essentially at 
all redshifts; however, at z = 6 (z = 10) the contribution of Pop III stars to the total 
galactic luminosity is always less than 5% for Mjjv < —17 {Muv < —16). The typical 
high-z galaxies closely resemble the GRB host galaxy population observed at lower 
redshifts, strongly encouraging the use of GRBs to detect the first galaxies. 
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The search for the most distant galaxies, located at the be- 
ginning of the cosmic dawn, is now entering its maturity. 
The last few years have witnessed a tremendous increase in 
the data available, and the number of candidates at red- 
shifts as high as z — 10, corresponding to only half a billion 
years after the Big Bang. This has been made possible by a 
combination of new technologies and refined selection meth- 
ods. In the first class of triggers, it is easy to acknowledge 
the role of the Hubble Space Telescope (HST). Thanks to 
dedicated surveys including the Hubble Ultra Deep Field 
and its predecessor, the Hubble Deep Field, we have been 
able to collect information on the luminosities and num- 
ber counts of galaxies located at the end of the reionization 
epoch. Immediately after, follow-up experiments performed 
with the newly installed Wide Field Camera (WFC3) , yield- 
ing sky images in the F105W (Y-band, 1.05/xm), F125W (J- 
band, 1.25/mi) and F160W bands (H-band, 1.60^im), have 
allowed to push the exploration to very faint (e.g. AB mag 



= 28.8 in the above bands) galaxies as remote as z — 10. 
In addition, the WFC3 crafted filters have considerably al- 
leviated the contamination problem due to interlopers and 
provided more precise photometric redshift estimates. The 
standard selection method applied to these survey data sets 
is based on the dropout technique introduced by Steidel 
et al. (1996) and later constantly refined and improved by 
several authors (e.g. Giavalisco et al. 2004, Bouwens et al. 
2007). Though this method has proved to be very solid in 
identifying high-redshift sources, it has the drawback that 
the exact source redshift cannot be determined with com- 
plete confidence. This uncertainty can be partly overcome 
by also using the longer wavelength infrared data, such as 
that provided by the Spitzer satellite; by building a more 
complete Spectral Energy Distribution (SED), the stellar 
mass, age and redshift of a given source can be constrained 
further. Other complementary techniques to search for dis- 
tant (z > 5) galaxies are also widely used, among which the 
narrow-band spectroscopy (Malhotra et al. 2005; Shimasaku 
et al. 2006; Taniguchi et al. 2005; Kashikawa et al. 2006) 
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aimed at detecting the Lya line, carrying a large fraction of 
the bolometric luminosity, is definitely the most established 
one. Such narrow-band searches have yielded the record- 
holding most distant galaxy at z = 6.96 (lye et al. 2006). 
Finally, another series of experiments involve searching for 
remote galaxies behind foreground galaxy clusters acting as 
magnification lenses (Schaerer & Pello 2005; Richard et al. 
2008; Bradley et al. 2008). Although these searches result in 
deeper magnitudes, their interpretation is hampered by the 
lens modelling and by the extremely narrow field of views, 
rendering it difficult to keep cosmic variance under control. 
As a final remark, we note that the most distant, spectro- 
scopically confirmed, cosmic object is a Gamma Ray Burst 
(GRB090423 at z = 8.2, Salvaterra et al 2009b; Tanvir et al. 
2009) . Although not a galaxy, the presence of this indicator 
implies that star formation was already well under way at 
those early epochs, thus further encouraging deeper galaxy 
searches. In addition, the GRB can be seen as a signpost of 
the underlying galaxy which could possibly be detected in 
the future knowing its exact position. Such a finding would 
be of the utmost importance as GRBs are mostly associated 
with star forming dwarf galaxies (Savaglio, Glazebrook & Le 
Borgne 2009) which are now considered to be the dominant 
sources of (re)ionizing photons at high redshifts (Choudhury 
& Ferrara 2007; Choudhury, Ferrara & Gallerani 2008). 

What have we learned from this wealth of experimen- 
tal results ? The most solid piece of information that can 
be determined from the data appears to be the luminos- 
ity function (LF) and, less robustly, its evolution. It is 
useful to briefly recap the present observational situation 
marching towards increasing redshift. Bouwens et al (2007) 
present a comprehensive view of galaxy candidates from the 
UDF/ACS/GOODS fields using NICMOS in the redshift 
range 2 = 4-6. They identify 1416 (627) V-dropouts (i- 
dropouts) corresponding to z ~ 5 (z ~ 6) down to an 
absolute UV magnitude of Muv ~ —17 with a LF de- 
scribed by a Schechter function with characteristic luminos- 
ity and faint-end slope given by M uv = —20.64 ± 0.13 and 
a = -1.66 ± 0.09 (M uv = -20.24 ± 0.19, a = -1.74 ± 0.16) 
respectively. The same group (Bouwens et al. 2008) has ex- 
tended the data analysis to include z ~ 7 z-dropouts (8 can- 
didates at z = 7.3) and J-dropouts (no candidates at z « 9). 
More recently, the installation of WFC3 on board the HST 
has triggered a new series of searches. Oesch et al. (2010) 
used data collected during the first-epoch WFC3/IR pro- 
gram (60 orbits) in the Y, J, H bands reaching a magnitude 
limit of AB?« 29(5ct). They identify 16 z-dropouts in the 
redshift range z = 6.5 — 7.5 from which they obtained a LF 
with (M uv = -19.91 ± 0.09, a = -1.77 ± 0.20), essentially 
confirming the previous findings while extending it to fainter 
luminosities (Muv ~ —18). Bouwens et al. (2010a) pushed 
the investigation to z = 8.0 — 8.5 by using 5 Y-dropouts. 
Finally, Bouwens et al. (2009) were able to identify three J- 
dropouts. If confirmed, these sources would be the most dis- 
tant objects detected so far. Similar studies using the same 
data has been performed by Bunker et al. (2010), who find 
a comparable number of z— and Y— dropouts. McLure et al. 
(2010) did not apply specific color cuts as in the previous 
works, thus finding a larger number of candidates; however, 
they pointed out that about 75% of the candidates at z > 6.3 
(100% at z > 7.5) allow a, z < 2 interloper solution. A re- 
cent analysis of the three HUDF and of the deep (~ 27.5 



AB mag), wide-area (~ 40 arcmin 2 ) WFC3 Early Release 
Science reveals 66 and 47 candidate galaxies at z = 7 and 
z = 8, respectively (Bouwens et al. 2010b). After carefully 
modelling the selection volume of each field and of the pos- 
sible contamination by spurious sources, the LF obtained 
from these data, while consistent with previous derivations 
of M uv and of the normalization, shows a steeper faint-end 
slope with a = -1.94 ± 0.24 and a = -2.00 ± 0.33 at z = 7 
and 2 = 8, respectively. 

Besides the LF, tentative information on the physical 
properties of these sources can be extracted from their SED, 
exploiting available Spitzer data (Eyles et al. 2005; Yan et 
al. 2006; Stark et al. 2009). In a recent study Labbe et al 
(2010b), based on follow-up Spitzer /IRAC observations, an- 
alyzed the SED of 12 z-dropout and 4 Y-dropout candidates. 
None of them is detected in the Spitzer/IRAC 3.6 /tm band 
to a magnitude limit of AB=26.9 (2a), but a stacking analy- 
sis reveals a robust detection for the z-dropout sample and a 
strong upper limit for the Y-dropout one. The stacked SEDs 
are consistent with a stellar mass of about 1O 9 M0, no dust 
reddening, sub-solar metallicity, and best-fit ages of about 
300 Myr, implying a formation epoch z £s 10. These results 
for the stacked sample should be compared with those ob- 
tained by Finkelstein et al (2010) who performed an object- 
by-object analysis and found similar ages but with a consid- 
erable spread, allowing ages as low as a few Myr. 

One of the major triggers to look for very high-z galax- 
ies is the quest for the reionization sources. The ionizing 
photon budget provided by the candidate high-z galaxies 
is often estimated by extrapolating their LF to lower lumi- 
nosities, a step that introduces a considerable uncertainty 
in the final determination. Having this in mind, it is still in- 
teresting to note that most studies tend to agree on the fact 
that the integrated UV specific luminosity for the detected 
galaxies at z = 7 — 8 falls short of accounting for the ion- 
izing power required to reionize the intergalactic medium. 
Of course, this conclusion is subject to at least two major 
unknown factors, these being the gas clumping factor (af- 
fecting its ability to recombine), and the escape fraction of 
ionizing photons (affected by dust and neutral hydrogen ab- 
sorption within galaxies). Additionally, the effect of poorly 
constrained ages and metallicities (including the presence 
of metal- free, massive Pop III stars), further complicate the 
calculation. 

In spite of the large experimental effort, surprisingly lit- 
tle attention has been devoted by modelers to the very high 
redshift universe. Most of the work has so far concentrated 
on a semi-analytical approach (Stiavelli, Fall & Panagia 
2004; Schneider et al. 2006; Bolton & Haehnelt 2007; Mao et 
al. 2007; Samui, Subramanian, Srianand 2009; Trenti et al. 
2010) to compute the luminosity function, number counts 
and emissivity evolution of high-z galaxies. Albeit quite fast 
and versatile, these methods cannot provide detailed infor- 
mation on the properties of the galaxies, often being based 
on simplified assumptions. Numerical dedicated simulations 
attempting to model galaxy populations beyond z = 5 — 6 
are also very scarce, with the partial exceptions constituted 
by the works by Nagamine et al. (2006) and Finlator, Dave 
& Oppenheimer (2007). 

Our approach is novel and different in spirit from all 
the previous theoretical ones. As our main aim is to model 
very high redshift reionization sources, we can afford smaller 
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simulation boxes, thereby reaching the high resolutions re- 
quired to resolve the dominant reionization sources - dwarf 
galaxies. Most importantly, though, we have implemented a 
careful treatment of metal enrichment and of the transition 
from Pop III to Pop II stars, along with a careful modelling 
of supernova feedback. Here we are interested in deriving 
the LF plus other observables from the simulations and to 
cast them in a form that can be compared directly with the 
available data or used to make new predictions for the James 
Webb Space Telescope (JWST). 



2 NUMERICAL SIMULATIONS 

For the present study we have performed a set of cosmo- 
logica j simulations using the publicly available code GAD- 
GET] 2 ] (Springel 2005) with an improved treatment of chem- 
ical enrichment as detailed in Tornatore, Ferrara & Schnei- 
der (2007, TFS07). A unique feature of the computation 
concerns the Initial Mass Function (IMF) of stars, which 
is taken to be different for Pop III and Pop II stars and 
depends on gas metallicity, Z. In brief, if Z < Z CI , the 
adopted IMF is a Salpeter law with lower (upper) limit 
of IOOMq (5OOM0) and only pair-instability supernovae 
(14OM < M < 260M Q ) contribute to metal enrichment 
(Heger & Woosley 2002). If Z > Z cr , the above limits are 
shifted to 0.1M Q (100M©), respectively; stars above 4OM 
end their lives as black holes swallowing their metals. Ac- 
cording to the canonical choice, we fix Z CT = 10~ 4 Z e . These 
two populations, to which we will refer to as Pop III and 
Pop II stars respectively, differ also for their metal yield and 
explosion energy. Complete details of the simulation can be 
found in TFS07. The simulation follows the production and 
transport of six different metal species, namely: C, O, Mg, 
Si, S, Fe, based on which the locally appropriate IMF is se- 
lected. The simulated volume has a linear (comoving) size 
L = lOh' 1 Mpc with N p = 2 x 256 3 (dark+baryonic) par- 
ticles, corresponding to a dark matter (baryonic) particle 
mass of M p = 3.62 x 1O 6 /i _1 M (6.83 x 10 5 /i _1 M©); the 
corresponding force resolution is 2 kpc. Our resolution does 
not allow us to track the formation of mini-halos, whose 
stellar contribution remains very uncertain due to radiative 
feedback effects (Haiman & Bryan 2006; Susa & Umemura 
2006; Ahn & Shapiro 2007, Salvadori & Ferrara 2009). The 
computation is initialized at z = 99 and carried on until 
2 = 2.5. Star formation is treated as a stochastic process, 
following Springel & Hernquist (2003) Basically, at a given 
time the star formation rate of a (multi-phase) gas particle 
is computed using a Schmidt-type law m = xm/t*, where 
xm is the mass of cold clouds providing the reservoir for 
star formation. Within the effective star formation model by 
Springel & Hernquist (2003), the star formation time-scale, 
is computed as t*(p) = ^(n/n,)" 1 / 2 , where t* = 1.5 Gyr 
and the density threshold for star formation, n* =0.1 cm -3 , 
are free parameters of the model chosen so as to reproduce 
the observed Kennicut relation within uncertainties. 



Throughout the paper, we adopt a ACDM cosmological model 
with parameters U M = 0.26, U A = 0.74, h = 0.73, U b = 0.041, 
n = 1 and trg = 0.8, in agreement with the 3-yr WMAP results 
(Spergel et al. 2007). 

2 www.mpa-garching.mpg.de/galform/gadget/ 



Supernova winds are treated as in the original model 
by Springel & Hernquist (2003); however, for simplicity and 
because the mass load and kinetic energy fraction are un- 
known parameters, we have given the winds from both popu- 
lations an initial velocity v w — 500 km s~\ which appears to 
be consistent with that derived from observations of high-z 
starburst galaxies. Wind particles are temporarily hydrody- 
namically decoupled until either (i) they have moved by a 
traveling length A = 2 kpc, or (ii) their density has decreased 
below 0.5 times the star formation density threshold. The 
metals are donated by star particles to the surrounding gas 
ones using a SPH kernel as described in TFS07. 

The gas photo-ionization and heating rates are calcu- 
lated at equilibrium with a background ionizing radiation 
due to the combined contribution of galaxies and quasars, 
taken from Haardt & Madau (1996), shifted so that the in- 
tensity at 1 Ryd is J v = 0.3 x 10~ 21 erg s -1 Hz -1 , in 
agreement with Bolton et al. (2005). Gas cooling and wind 
treatment details are the same as in TFS07. 



3 LUMINOSITY FUNCTION 

As a first check of the simulation results we will compute the 
evolution of the LF of the simulated galaxies and compare 
it with available data. To this end we make a number of 
physical assumptions that are discussed in the following. 

For each galaxy at a redshift z, the luminosity at rest- 
frame wavelength, A, is the sum of the contribution of its 
Pop II and Pop III stars. The SED of Pop II stars is com- 
puted by running the Starburst99 code (Leitherer et al. 
1999, Vazquez & Leitherer 2005) using the metallicities and 
stellar ages appropriate for the galaxy under consideration 
as obtained from the cosmological simulation. The SEDs of 
massive (metal-free) Pop III stars have been computed by 
Schaerer (2002) including the effect of nebular emission. For 
both populations the IMF is taken according to the same 
prescription used in the simulation and described in Sec. [2] 
The total luminosity of a galaxy is then given by 



, rll , rill ,11/ II v-i.V 11 , jiiiiV 111 1 

L x = L x + L x = l x (r , Z)M t +l x M* t 



(1) 



where l\ l {r ll ,Z) is the SED template for Pop II stars with 
mean age r 11 and metallicity Z corresponding to a continu- 
ous star formation rate of 1 Mq yr _1 ; Af* is the Pop II star 
formation rate. The SED template of Pop III stars, whose 
star formation rate is M, , is l\ u ; r ln = 2.5 x 10 6 yr is the 
mean lifetime of massive Pop III stars (Schaerer 2002). We 
have implicitly assumed that the emission properties of Pop 
III stars are constant during the short lifetime of these mas- 
sive metal- free stars. Next, the luminosity is converted into 
absolute AB magnitudes, Muv, where the suffix UV refers 
to the wavelength A = (1600, 1350, 1500, 1700, 1500, 1600) 
A for z = (5,6,7,8,9,10), respectively. Though this wave- 
length choice is chosen to be consistent with the observa- 
tional data we compare with, a different value would not 
affect the results in any sensible way, given the flatness of 
the spectrum in this short wavelength range. 

The LF of galaxies at any redshift z is obtained by 
counting the galaxies with a given absolute magnitude in 
each magnitude bin and dividing the final result by the 
total volume of the simulation and bin size (0.5 mag). 
We perform this procedure for the following six redshifts 
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Figure 1. The UV LF of galaxies at the different redshifts shown in each panel. Observational data (and upper limits) from HUDF are 
taken from Bouwens et al. (2007) for z = 5,6; Oesch et al. (2010; squares), McLure et al. (2010; triangles) and Bouwens et al. (2010b; 
circles) for 2 = 7; Bouwens et al. (2010a; squares), McLure et al. (2010; triangles) and Bouwens et al (2010b; circles) for 2 = 8 and 
Bouwens et al. (2009) for z = 10; they are shown as circles (arrows). The histograms show the simulated LF with error bars representing 
Poisson errors. Dotted lines are the Schechter function fits to the LF; the vertical short (long)-dashed lines mark the sensitivity limit of 
JWST (HST/WFC3). 




Figure 2. Maps showing the distribution of galaxies with an observed flux larger than the JWST sensitivity limit of 1 njy in the J-band 
(left panel) and H-band (right). The maps are a 2D cut of a 3D image produced by stacking simulation snapshots between z ~ 7.6 — 11.6. 
The vertical color bar gives the redshift of the object, the size of the galaxies scales with their flux in the range 1 njy < T < 24.4 njy. 
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Figure 3. the color excess, E(B — V), as a function of the stellar 
mass for z 7 galaxies obtained from the simulation snapshot. 
The dark (light) region shows the color excess for 65% (95%) of 
the galaxies. 



z = (5, 6, 7, 8, 9, 10). The results are shown in Fig. [T] as solid 
histograms, where the error bars represent the Poisson er- 
ror on the number of galaxies in each magnitude bin. These 
theoretical LFs are then compared to the experimental ones 
collected from the various analyses of the HUDF. For z — 10 
we show the upper limits on the LF obtained from the three 
available candidates identified by Bouwens et al. (2009). 

Let us now analyze the results shown in Fig. [T] in more 
detail. It is clear that the luminosity range sampled by the 
observations and our predictions is only partially overlap- 
ping. This is because on one hand, even the exquisite sensi- 
tivity of WFC3 is not sufficient to properly sample the faint- 
end of the LF (Muv si — 18); on the other hand, our simula- 
tions, which are specifically designed to properly resolve the 
very first galactic units in a relatively small volume, lack the 
most massive, rare objects which comprise the bright end of 
the LF. In spite of these shortcomings, we consider it a re- 
warding success that the amplitudes of the theoretical and 
experimental LFs match almost perfectly, and at the same 
time, have quite similar slopes at all redshifts for which data 
are available. This is even more striking as no attempts have 
been made to fit or adjust the theoretical curves to the ob- 
served LF, i.e. they have been computed directly from the 
simulation output with no free extra parameters. 

Our results suggest two clear trends. First, the LFs shift 
towards fainter luminosities with increasing redshift, mim- 
icking a pure luminosity or density evolution. This is quite 
consistent with the trend of an increasing M'^ v with red- 
shift, found in the data by several groups (see the extended 
discussion in Ouchi et al. 2009) , preferring a pure luminosity 
evolution. Second, the faint-end slope of the simulated LF 
does not vary (within errors) from z — 5 to z = 10, main- 
taining an almost constant value of a = —2.0. This value is 



slightly larger than the one derived from the data (see Sec. 
[TJ at z — 5 — 6 but perfectly consistent with the recent de- 
termination of a for the z = 7 — 8 LF obtained by Bouwens 
et al. (2010b). This behavior is most likely produced by a 
combination of the halo mass function evolution, feedback 
effects and evolving stellar populations. 

Clearly the faint-end slope will be better constrained by 
forthcoming facilities such as the JWST; the above success- 
ful test of our model allows us to make reliable predictions 
for surveys that will be performed with such instruments. 
For a deep exposure of 10 6 s, JWST is expected to reach 
a photometric sensitivity of about 1 njy (lOcr) allowing an 
investigation of the faint-end of the LF predicted by our 
simulations at least up to z — 9 (see Fig. [T]). This will be 
particularly exciting because it will very likely unveil the 
physical properties of these objects which are now thought 
to be the main reionization sources, as we will discuss in 
detail in Sec. [5] We want to stress here that our simula- 
tion allow us to resolve properly (i.e. with more than 200 
DM particles) the galaxy population that will be observed 
by JWST, being their total halo masses > 10 9 M©. This, 
toghether with the carefull threatment of the chemical feed- 
back, makes our predictions particularly robust. 

To give a visual idea of the high-z (z > 7.5) galaxy 
population in a typical JWST field (of size 5.3 arcmin) we 
have produced the maps shown in Fig. [2]which demonstrate 
the wealth of objects that JWST will be able to detect up 
to z > 11. The maps are a 2D cut of a 3D image produced 
by stacking all the available simulation snapshots between 
z ~ 7.6 — 11.6, and filling in the gaps between snapshots 
by randomly rotating the orientation axes in all the three- 
dimensions simultaneously. 

As it has emerged from the overview of the observa- 
tional results given in the Introduction, dust seems to have 
only a minor effect on shaping the SED of faint high-redshift 
galaxies. This is easily understood, as we will see later on, as 
a result of the relatively low stellar masses and metallicities 
characterizing these objects. Dust can however become more 
important in luminous Lyman Break Galaxies (LBGs) and 
Lyman Alpha Emitters (LAEs) as many studies have shown 
(Lai et al. 2007; Atek et al. 2008; Nagamine, Zhang & Hern- 
quist 2008; Dayal et al. 2008, 2009, 2010, 2011; Finkelstein 
et al. 2009). 

The dust enrichment of galaxies can be followed self- 
consistently by post-processing our simulation through the 
model introduced by Dayal et al. (2010), using which, we can 
quantitatively check the previous statement. In brief, the 
model assumes that dust is produced by supernovae (we ne- 
glect here the contribution of AGB stars, though see Valiante 
et al. (2009) who found a significant contribution from these 
objects for population ages > 150 — 200 Myr) and taking into 
account three processes: (a) dust forms in the expanding 
ejecta with a yield per SNII of 0.54M© (Todini & Ferrara 
2001; Nozawa et al. 2003; Bianchi & Schneider 2007), (b) 
SNII destroy dust in the ISM they shock to velocities > 100 
km s _1 , with an efficiency of 0.12 (McKee 1989), and (c) 
a homogeneous mixture of gas and dust is assimilated into 
star formation (astration) . Once the dust mass is calculated 
for each galaxy in the simulation, we translate this into a 
value of E(B — V) using the appropriate SN dust extinction 
curve given by Bianchi & Schneider (2007) as explained in 
Dayal et al. (2010). The resulting values of E(B - V) for 
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galaxies at z ~ 7 are shown in Fig. [3] as a function of the 
stellar mass. Many galaxies, especially the smallest ones, are 
almost dust-free, and none of them shows a dust reddening 
value larger than E(B — V) — 0.009. This evidence allows 
us to safely neglect the effects of dust on the UV LF. 



4 PROPERTIES OF FIRST GALAXIES 

In addition to predicting the global evolution of the LF, a 
major strength of our study is that it makes possible to ex- 
tract the physical properties of the high-z galaxies which are 
a part of the faint-end of the LF. We start by concentrating 
on the stellar ages, reported in Fig. [4] as a function of the 
UV magnitude, from which the following main conclusions 
can be drawn. On average, fainter (and less massive) ob- 
jects tend to be younger at all redshifts, with typical ages 
for observable objects in the range 200-300 Myr at z — 5 
and 80-130 Myr at z — 7 — 8, with the caveat that a consid- 
erable age spread exists at all luminosities. These ages imply 
that these galaxies started to form stars as early as z = 9.4, 
clearly suggesting that their UV light might have influenced 
the cosmic reionization history. From a different perspective, 
we note that stellar ages of ~ 100 — 150 Myr are expected 
for galaxies at the sensitivity limit of WFC3 for z = 7 — 8, 
in agreement with the observational estimates obtained by 
Labbe et al. (2010a) and Finkelstein et al. (2010). 

A second important key physical parameter of pristine 
galaxies is their stellar and gas metallicity. In the follow- 
ing we will only discuss the stellar metallicity, keeping in 
mind that the two closely match each other. Somewhat sur- 
prisingly (but not unexpectedly) even the faintest galaxies 
appear to be already enriched to remarkable levels: at the 
JWST sensitivity threshold, we find Z > 0.03Z© at all red- 
shifts; galaxies identified by HST (Muv < —18) system- 
atically show metallicities in excess of Z — 1/IOZq even 
at z — 7 — 8. Thus we come to the interesting conclu- 
sion that even at these early epochs, the self-enrichment, 
due to the metals produced, following the first star forma- 
tion episodes is able to increase the metal abundances of 
such small objects to levels comparable to their present-day 
counterparts (e.g. the Magellanic Clouds). In addition, such 
high mean metallicites could in principle preclude the for- 
mation of Pop III stars according to the critical metallicity 
scenario (Schneider et al 2002, 2003; Schneider & Omukai 
2010) which predicts Z cr — 10~ 5±1 Zq as the upper limit for 
the formation of Pop III stars. However, as we will discuss in 
more detail later, inside these early structures, small pocket 
of (quasi) pristine gas may survive in which a relatively tiny 
amount of Pop III stars continue to form as pointed out by 
TFS07 and Jimenez & Haiman (2006). In brief, the scenario 
leads to the concept of a "PopIII wave", i.e. the physical 
phenomenon by which in each galaxy the formation of stars 
below the critical metallicity is progressively segregated to- 
wards the external regions of the galaxy, where almost un- 
polluted regions are still present. Until this process comes to 
an end (when metal pockets produced around the first stars 
forming regions reach a considerable volume filling factor in 
the system), PopIII and PopII formation modes coexist in 
the same galaxy. One might naively think that this process 
can take longer in large galaxies than in small ones, but it 
is (roughly) true only to zero-th order approximation, as in 



lower mass galaxies the efficiency of star formation is also 
depressed and the amount of time required, once scaled by 
the baryonic mass of the system, is not dramatically differ- 
ent from that of larger systems. 

Taken together, Fig. [4] and [5] provide a first guess of 
the properties of a typical high-z galaxy with given absolute 
UV magnitude Muv- This can be useful when comparing 
the observed photometric data with synthetic SEDs. 

Fig. [5] shows the relation between the SFR and stellar 
mass of galaxies at z = 7. Simulated objects follow an almost 
linear relation with significant spread, i.e. an almost con- 
stant specific star formation rate. This trend closely matches 
the one found in the analysis of stacked SEDs of WFC3 z- 
dropout by Labbe et al. (2010b), although a single object 
analysis reveals large errors in the determination of both 
the stellar mass and the SFR for these objects (Gonzalez et 
al. 2010). We note that simulated galaxies tend to have a 
slightly higher SFR (for a fixed stellar mass) than expected 
by the extrapolation of the relation to smaller objects. How- 
ever, in a re-analysis of the observed z = 7 sample including 
the effect of nebular continuum and line emission along with 
that of dust absorption, Schaerer & de Barros (2010) find 
smaller stellar masses and larger SFR with respect to previ- 
ous works. Our simulated galaxies lie somewhat in between 
these two observational estimates. 

Fig. [7] shows the evolution of the galaxy properties as a 
function of redshift for galaxies with different stellar masses. 
The mean age of stellar population for all galaxies in the 
simulation is found to decrease with redshift oc (1 + z)~ 2 ; 
however, the age spread increases for less massive galaxies at 
any redshift. The average stellar metallicity is also growing 
with time for M* > 10 7 A#q, but show a much flatter (almost 
constant) evolution in the smallest objects, leveling at about 
1/25 Zq for the tiniest star-forming galaxies. However, even 
for dwarf galaxies with M* « 10 5 M© one can find individual 
objects enriched up to 1/10 of solar metallicity already at 
z = 10. This mass-dependent metallicity evolution is proba- 
bly caused by the different ability to retain metals deposited 
by supernova explosions of the massive and dwarf popula- 
tions. As the potential wells of the latter one are shallower, 
metals escape easily into the intergalactic medium thus set- 
ting an upper limit to the amount of metals than can be 
kept in their main body (Mac Low & Ferrara 1999). An 
obvious implication is that IGM metals preferentially come 
from small and common objects, thus resulting in a more ho- 
mogeneous enrichment (Ferrara, Pettini & Shchekinov 2000; 
Madau, Ferrara & Rees 2001). 



5 REIONIZATION SOURCES 

Having obtained the luminosity function and characterized 
some of the physical properties of the current high-z candi- 
date galaxies, we now turn to an analysis of their role in the 
reionization process. 

The rate of ionizing photons from the j-th galaxy is 
given by 

ritonj = fesc [Q u Mf + Q m Mf l r m ] (2) 

where Q 11 (Q™) is the ionizing photon flux for Pop II (Pop 
III) stars, which of course is a function of the IMF, metal- 
licity and stellar age of any given galaxy. Finally f esc is 
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Figure 4. Mean stellar age of galaxies located at different redshifts (see labels) as a function of their absolute UV magnitude. The dark 
(light) shaded area show 68% (95%) of the galaxies in the magnitude bin. The vertical short (long)-dashed lines mark the sensitivity 
limit of JWST (HST/WFC3). 
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Figure 5. Mean metallicity of galaxies located at different redshifts (see labels) as a function of their absolute UV magnitude. The dark 
(light) shaded area show 68% (95%) of the galaxies in the magnitude bin. The vertical short (long)-dashed lines mark the sensitivity 
limit of JWST (HST/WFC3). The filled circle refer to the GRB 050904 host galaxy at z = 6.3. 
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Figure 6. Stellar mass-SFR relation for simulated galaxies at 
z = 7 (black points) compared with recent estimates for ob- 
served z-dropout as derived by different authors: triangles are 
from Schaerer &; de Barros (2010), open squares from Gonzales 
et al. (2009), and filled squares for the mean values derived by 
Labbe et al. (2010b). The dotted line shows the empirical relation 
between the two quantities obtained by Labbe et al. (2010a). 



the escape fraction of ionizing photons. For simplicity, we 
assume a redshift-independent value of / esc for the two stel- 
lar populations. The ionization rate provided by galaxies in 
the simulation box per unit comoving volume at redshift z, 
Nion(z), is then given by the sum over all galaxies at that 
redshift divided by the volume of the simulation. 

The actual ionization rate must also include galaxies 
that are too rare to be caught in our relatively small sim- 
ulation volume. To account for this correction, we add the 
ionization rate due to bright galaxies by integrating the ob- 
served LF (or the upper limits) from the luminosity of the 
brightest galaxy in the output down to very low magnitudes 
(Muv = —25). The ionizing photon flux is obtained using a 
SED derived assuming that rare galaxies have the same stel- 
lar age and metallicity as the brightest simulated one. The 
extra contribution of the unaccounted bright-end of the LF 
is found to be at most 10% of the total ionizing photon 
emission. 

The evolution of the total specific ionization rate, Ni OTt , 
is plotted in the top panel of Fig. [5] (/esc = 0.2) along with 
the same rate due to galaxies detected in the HUDF, N^ n , 
or detectable by JWST, N/ on . The ratios between (N^ n , 
Nf on ) and Ni 0n are also shown for clarity in the bottom 
panel of the same Figure. HST is now resolving the sources 
that provide ~ 1/3 of the ionizing photon budget at z = 5 
and ~ 20% at z = 7 — 7.5. This results is consistent with 
the limits set by the observed very steep faint-end slope of 
the z = 7 — 8 LF (Bouwens et al. 2010b) and with the 
estimate of recent semi-analytical models (Choudhury et al. 
2008, Trenti et al. 2010). 

At the sensitivity limit of JWST, it will be possible to 
detect the bulk of ionizing sources up to z ~ 7.3, but at 
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higher redshifts most of the ionizing photons will still be 
produced by sources that are too faint to be detected even 
by JWST. 

The total ionization rate density N ion (z) should then be 
compared with the recombination rate density of the IGM, 
Nrec(z), given by (e.g. Madau, Haardt & Rees 1999) 



l + z 
7 



_1 Mpc" 



(3) 



where (nji) is the mean comoving hydrogen density in the 
Universe and (t rec ) is the volume- averaged recombination 
time for ionized hydrogen with an effective Hn clumping fac- 
tor Chii = {n'j II i) I ' (nun) 2 ■ The recombination rate density 
is shown in the top panel of Fig. [5] with dotted lines for dif- 
ferent value of the clumping factor Chii- For f esc — 0.2 and 
Chii = 10, the balance between ionization and recombina- 
tion is obtained at z ~ 6.8. For f e3C = 0.1 Ni OU = Nrec at 
2 = 6 assuming Chii = 10. 



6 LIGHT FROM POP III STARS 

Another piece of useful information than can be extracted 
from the simulation outputs is the relative fraction of nor- 
mal (Pop II) and massive, metal-free (Pop III) stars. There 
are several questions to which we can provide quantitative 
answer: (i) do some of the current candidates contain Pop III 
stars ? (ii) in that case, what fraction of their UV luminosity 
is powered by them ? (iii) how is this fraction dependent on 
their Muv luminosity ? (iv) is there a clear observational 
signature imprinted by Pop III stars ? 

The answer to the first question is straightforward: 
having analyzed the stellar populations of the simulated 
galaxies present at four observationally relevant redshifts, 
2 = (5, 6, 8, 10), we find that a fraction 0.07-0.19 (depending 
on z) of the galaxies contain at least some Pop III stars. We 
should not emphasize too much on the exact values of this 
Pop Ill/Pop II galaxy ratio as fluctuations in galaxy mass 
and star formation rate might introduce a very large disper- 
sion. The most robust physical quantity to understand the 
relative importance of the two populations is the ratio of 
Pop III-to-Pop II star formation rates which is a decreasing 
function of time (see Fig. 1 of TFS07) never exceeding 10 -3 
below z — 10. 

More relevant is question (ii) above, whose answer 
can be obtained by inspecting Fig. [9] showing the ratio, 
Luv,iii / Luv of UV luminosities contributed by Pop III and 
Pop II stars as a function of Muv- Light from Pop III stars 
becomes progressively more important towards fainter ob- 
jects. This confirms previous findings (e.g. Schneider et al. 
2006) that Pop III stars preferentially form in low-cr peaks 
rather than in larger galaxies whose gas has already been 
enriched by several stellar generations. At the same time, 
we do not find galaxies containing only Pop III stars: this 
occurrence is made unlikely by the short lifetimes of such 
very massive stars. At z = 6 (z = 10) the contribution of 
Pop III stars to the total luminosity is always less than 5% 
for M uv < -17 {Muv < -16). 

Among the candidate galaxies detected so far, we find 
that Pop III stars contribute less than a few percent to 
the total galaxy luminosity. Even at the detection limit 
of JWST, no Pop Ill-dominated galaxies (> 50% of total 
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Figure 8. Upper panel: Redshift evolution of the total specific 
ionization rate (solid line). The short (long)-dashcd line corre- 
sponds to galaxies detectable by JWST(HST/WFC3); dotted 
lines show the specific IGM recombination rate for different val- 
ues of the clumping factor Chii- Bottom panel: Fraction of ioniz- 
ing photons coming from galaxies identified by JWST and in the 
HUDF. An escape fraction f esc = 0.2 has been assumed. 



luminosity) will be found, due to their extreme faintness 



> -14. at z 



6). However, some objects having 



~ 10% of their luminosity powered by Pop Ills are present 
at z — 6 at the 1 njy sensitivity level reachable by the 
JWST. 

Recently, many authors (Shapley et al. 2003, Nagao et 
al. 2008, di Serego Alighieri et al. 2008) have searched for the 
Hell 1640A emission line in the spectra of the so-called dual 
emitters, i.e. high-2 galaxies showing strong emission in both 
Lya and Hell lines. The Hell line is usually taken as a well- 
defined signature of Pop III stars. Up to now, these searches 
have given negative results. Motivated by these attempts, we 
have computed the rest-frame equivalent width of the Hell 
line for the objects in our simulation box as £W(HeII) = 
L(HeII)/L ^ where L(HeII) = 5.67 x 10 _12 Q(He + ) is the 
luminosity of the Hell line in erg s _1 , Q(He + ) is the rate of 
He ionizing photons, and L ° is the continuum luminos- 

b 1 ' 1460A 

ity at A = 1460A in erg s A .In particular, we focus 
here on z — 6 which might be more easily accessible to 
present or future observations. Regretfully, the perspectives 
of direct detection of PopIII stars through this technique do 
not appear as very promising. For objects detectable in the 
JWST (HUDF) deep field survey, the expected Hell rest- 
frame equivalent width is < 0.5 A (< 0.1 A). Such small 
EWs will be very difficult to detect. The EW increases to 
more accessible values of about 10 A only if much fainter 
objects (Muv = — 13) could be observed. We have to un- 
derline that the above discussion implies that we cannot set 
limits on the total cosmic SFR in Pop III stars using the 
results of dual emitters searches given that the bulk of Pop 
III stars may be 'hidden' in galaxies much fainter than any 
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Figure 9. Fraction of the total luminosity Ltjv due to Pop III 
stars Li/y jfj as a function of the absolute UV magnitude Muv- 
The vertical short(long)-dashed lines mark the sensitivity limit of 
JWST (HST/WFC3). Absolute UV magnitudes are computed as 
in Fig.[T] 



present (and probably future) survey can detect (Scanna- 
pieco, Schneider & Ferrara 2003). 

In conclusion, Pop III stars are found to form essen- 
tially at any redshift and in < 20% of the galaxies. However, 
their contribution to the total galaxy luminosity is very low, 
apart from the very faint objects and their detection could 
be extremely difficult even with the next generation of space 
telescopes. 



7 EARLY GAMMA-RAY BURSTS 

Long Gamma-Ray Burst (GRBs) are powerful flashes of j- 
rays that are observed with a frequency of about one per 
day over the whole sky. The 7— ray emission is accompanied 
by a long-lasting tail, called afterglow, usually detected over 
the whole electromagnetic spectrum. Their extreme bright- 
ness easily over-shines the luminosity of their host galaxy 
and makes them detectable up to extreme high redshifts, 
as shown by the discovery of GRB 090423 at z = 8.2 (Sal- 
vaterra et al. 2009b, Tanvir et al. 2009). Metal absorption 
lines can often be identified in their afterglow spectra, al- 
lowing a study of the metal (and dust) content of the en- 
vironment in which they blow. Finally, once the afterglow 
has faded, follow-up searches of the GRB host galaxy be- 
come possible. At low redshifts, GRBs are typically found 
in blue, low-metallicity dwarf galaxies with stellar masses 
M* ~ 10 8-9 Mq and high specific star formation rates 
(Savaglio et al. 2009). These objects closely resemble the 
properties of high-z galaxies identified in our simulations, 
whose mean specific star formation rates (SSFR) are w 8—10 
Gyr -1 , albeit associated with a large spread (for M* ^ 10 s , 
SSFR=1.5-7 Gyr" 1 at z = 6 - 8). This suggests that high-2 
GRBs can be used as signposts of the same faint galaxies 



that provide the bulk of the ionization photons in the early 
Universs[f| 

. Moreover, the study of their afterglows can provide 
new hints about the metal (and dust) content of the parent 
galaxy. Finally, it has been proposed that even Pop III stars 
may eventually blow as powerful GRBs (Fryer, Woosley & 
Heger 2001; Yoon, Langer & Norman 2006; Hirschi 2007; 
Komissarov & Barkov 2010). The observation of these Pop 
III-GRBs may provide a valuable way to detected the elu- 
sive, short living first stars. We will discuss these points in 
the following. 



7.1 GRBs as metallicity probes 

Fig. [5] shows the mean stellar metallicity of the simulated 
sample of high-z galaxies as a function of their absolute UV 
magnitude at different redshifts z ^ 5, i.e. the luminosity- 
metallicity relation for such objects. Although the experi- 
mental determination of such a relation would be of the ut- 
most importance in testing theoretical predictions, in prac- 
tice such effort is hampered by the extreme difficulty to infer 
Z from lines in the galaxy spectrum with available facilities. 

This problem can be considerably alleviated if a GRB 
could be found inside of one of these remote galaxies. In 
this case absorption features produced by heavy elements 
dispersed in the interstellar medium (ISM) surrounding the 
GRB would leave a characteristic and recognizable imprint 
on top of the GRB afterglow spectrum. Moreover, there is 
now a good agreement on that GRBs are usually hosted in 
relatively small, star- forming galaxies (Savaglio et al. 2009): 
these requirements would make the dwarf high- z galaxies 
which, according to our findings, are the dominant popula- 
tion during the first cosmic billion year, optimal candidate 
hosts. 

The only high-quality spectra burst at z > 6 available 
so far, GBR 050904 at z = 6.3 (Kawai et al. 2006), shows 
the expected metal absorption features, witnessing the pres- 
ence of metals at the same redshift of the GRBs (Kawai et 
al. 2006). From this data, and with the further assumption 
that the measured sulfur ([S/H] = — 1.3±0.3) is a good proxy 
for metallicity, one can determine the metal content of the 
host galaxy. Berger et al. (2007) attempted a detection of the 
GRB host galaxy with HST and Spttzer, and were able to 
set an upper limit to the host luminosity of Muv —20.3, 
i.e. Muv ^ M'uv at 2 = 6. By combining the HST and 
Spttzer upper limits, they set an upper limit to the stellar 
mass of the galaxy, M* < few x 10 9 Mq. This constraint is 
consistent with the predictions presented in Fig. [5] (notice 
the data point corresponding to the Berger et al. 2007 deter- 
mination) , which might indicate that at z fa 6 a galaxy with 
Z ~ 0.05^0 is on average 2-3 magnitudes fainter than the 
detection limit of the experiment; in addition, such metal 
abundance level is quite typical for galaxies located at that 
epoch. 



3 Note that GRBs may be biased tracers of the cosmic star forma- 
tion rate (e.g. Daigne et al. 2006; Salvaterra & Chincarini 2007; 
Salvaterra et al. 2009a). However, if this is related to the existence 
of a metallicity threshold for GRB formation of Z < 0.3 Zq as re- 
quired by the collapsar model (e.g. MacFadyen &; Woosley 1999) 
then we expect that high-z GRBs can be fair cosmic tracers. 
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GRB 090423 (Salvaterra et al. 2009b, Tanvir et al. 2009) 
has been recently detected at z = 8.2. The intrinsic proper- 
ties of the burst (both of the prompt and afterglow phases) 
are similar to those observed at low/intermediate redshifts, 
suggesting that the progenitor and the medium in which the 
burst occurred are not markedly different from those of low- z 
GRBs (Salvaterra et al. 2009b, Chandra et al. 2010). Chary 
et al. (2009) observed the field of GRB 090423 for 72 hours 
using Spitzer/TRAC at 3.6 /im, looking for its host galaxy. 
The observation was performed ~ 46 days after the GRB 
trigger, (corresponding to about 5 days in the burst rest- 
frame). A weak source was detected at the location of the 
GRB afterglow with Lab = 27.2 ± 0.3. The observed flux is 
consistent with the power-law decay of the GRB afterglow, 
suggesting that the source is still contaminated by the GRB 
emission and implying a limit on the absolute magnitude 
of the host of A^ 3goo ; > —19.96. From Fig. 3, we expect 
that the host galaxy of GRB 090423 should be enriched at a 
level of a few percent solar. This relatively high metallicity 
may explain the high equivalent neutral hydrogen column 
density measured in the X-ray afterglow, although a wide 
range of systematic effects or the contamination by low-z 
intervening absorption systems may be a more valuable, al- 
ternative explanation for the observed absorption (Chandra 
et al. 2010). 

These examples nicely show how GRBs represent an 
unique tool to study the high-z Universe. In particular, 
the observation of their optical-NIR afterglow may allow to 
eventually study the evolution of the mass-metallicity rela- 
tion up to very high-z. To achieve this goal, high resolu- 
tion and good signal-to-noise afterglow spectra are required 
as soon as possible after the GRB detection in 7-rays. As 
nicely demonstrated by GRB 090423, GRBs are easily de- 
tectable well beyond any other astrophysical object. Indeed, 
up to ~ 5% of all GRBs detected by the Swift satellite are 
expected to be at z > 6 (Salvaterra &: Chincarini 2007, 
Salvaterra et al. 2009b). Future missions (e.g. EXIST, XE- 
NIA, SVOM) will rapidly increase the high-z GRB sample 
(Salvaterra et al. 2008) allowing a statistical study of GRB 
hosts in the high-z Universe and a direct check of the galaxy 
metal enrichment history at those early epochs. In particu- 
lar, EXIST, thanks to its 1.1m optical-NIR telescope, will 
be able to take the GRB afterglow spectrum only 300 s after 
the trigger, allowing an on-board direct measure of the red- 
shift and the identification of metal absorption lines when 
the afterglow is still sufficiently bright even for high-z bursts 
(Grindlay et al. 2009). 

7.2 GRBs from Pop III stars 

It has been widely discussed (see Bromm & Loeb 2007 for a 
review) whether GRBs can arise from the collapse of mas- 
sive, metal-free stars. While the large envelopes of these ob- 
jects may suppress the emergence of relativistic jets out of 
their surface, the weak winds expected for low-metallicity 
stars can prevent angular momentum loss during their evo- 
lution, producing the rapidly rotating central configurations 
needed to produce a GRB. In the following we will assume 
that massive Pop III stars produce GRBs. 

Excluding Pop III progenitors in the mass range 
[140,260] M Q , that will explode as Pair Instability Super- 
novae (PISN, Heger & Woosley 2002) leaving no remnant, 



we have two possible channels through which GRBs can oc- 
cur: (i) 100-140 M (Yoon et al. 2006, Hirschi 2007) and 
(ii) 260-500 M (Fryer et al. 2001, Komissarov & Barkov 
2010). In the latter case, since both the luminosity and the 
duration are thought to be proportional to the black hole 
mass, an extremely bright and long GRB is expected (Fryer 
et al. 2001), while more typical luminosities and durations 
may be expected in the case of smaller progenitors (Hirschi 
2007). 

We can estimate a strong upper limit for the rate of Pop 
III-GRB detections as follows. Suppose a fraction foRB of 
all Pop III stars with masses in the ranges discussed above 
produce GRBs with a typical beaming angle 6, and that all 
Pop III-GRBs are detectable by present-day satellites given 
the their extreme brightness. Then the observed rate of Pop 
III-GRB, Rgrb(> z), is given by 

Rgrb{>z) = 5.5xl0- 3 (^) 2 f dzfaRBV^j 1 ^) 

where p III (z) is the total Pop III star formation rate at 
redshift z, r] = J^' m " (f)(m,)dm,/ m*0(m*)dm* is 
the number of Pop III stars with masses in the range 
[M m i n , Mmax] given the assumed IMF <f>, and V is the cos- 
mological volume per unit solid angle. The factor (1 + z) _1 
accounts for time dilation effects due to redshift. The com- 
puted rate of Pop III-GRB at z > 6 is then Rqrb ~ 
7.5(2.5)/grb(6>/6°) 2 yr" 1 sr" 1 for GRB progenitors in the 
mass range 100-140 M e (260-500 M ). This rate is of the 
same order of that expected for normal (i.e. Pop II/PopI pro- 
genitors) GRBs (Salvaterra & Chincarini 2007, Salvaterra 
et al. 2009b). However, we stress here that very likely 
fGRB <C 1, resulting in lower detection rates. The detec- 
tion of one of such Pop III-GRBs might represent the most 
promising way to directly detect the very first stars to have 
formed. 



8 CRITICAL DISCUSSION 

The present results allow to build a coherent and quantita- 
tive description of the properties of elusive high-redshift, 
possibly primordial, galaxies. They are also very useful 
to interpret the data coming from deep surveys as the 
HST/WFC3 and future ones. However, there is considerable 
room for improvement left by our study. In the following we 
would like to elaborate on the uncertainties and shortcom- 
ings of our findings. 

We first note that resolving the dwarf galaxy popula- 
tion and following the PopIII transition process along with 
the large variety of physical processes implemented in the 
simulation limits the size of the cosmic volume that can be 
simulated. Resolution is certainly an important issue, as it is 
well known (see, e.g. Governato et al. 2010) to affect the sim- 
ulated star formation rates and cause the loss of sub-galactic 
structures. The dependence of the results from resolution 
has been presented an analyzed for the same set-up of the 
present simulations in Tornatore et al. (2007) to which we 
refer the interested reader (see Fig. 1 of that paper). Res- 
olution might also alter the details of the "Pop III wave" 
evolution, since galactic substructure allows star formation 
to occur at the edges of galaxies as well. The PopIII-PopII 
transition is also very dependent on the assumed IMF of 
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PopIII stars. Our conclusions are valid under the assump- 
tion that PopIII stars were very massive (M ^ 100 M©) and 
the first metal production is driven by the explosion of pair- 
instability SNe (see Schneider et al. 2006 for alternatives). 

Our box is also too small to properly describe cosmic 
reionization, let alone that we are not even attempting to 
properly treat radiative transfer. These issues have been al- 
ready addressed in previous works of our group; as already 
stated we are concerned here with the properties of high-z 
galaxies, which are presumably more affected by their inter- 
nal physics rather than by the environment, as we explain 
below. Comparing the simulated volumes with the observed 
ones is very challenging as considerable uncertainty exists on 
the latter (see discussion in Appendix B of Bouwens et al. 
2009). However, as it could be induced from the extension 
of the LF towards the most luminous and rarest objects at 
z=7, we estimate that the volume sampled by experiments 
should be about 30 times larger than our simulated one. 

The next caveat comes from the fact that we have ne- 
glected the effects of minihalos (virial temperature < 10 4 
K). Our resolution does not allow us to track the forma- 
tion of such objects, whose stellar contribution remains 
very uncertain due to radiative feedback effects (Haiman & 
Bryan 2006; Susa & Umemura 2006; Ahn & Shapiro 2007, 
Okamoto, Gao & Theuns 2008, Salvadori & Ferrara 2009). 
The presence of such small collapsed structures, if able to 
form stars, could alter the reionization history to some ex- 
tent and increase the number counts of high-redshift galax- 
ies, if detectable. As far as reionization is concerned, it has 
already been shown by Choudhury & Ferrara (2007) that 
acceptable fits to all relevant reionization data can be ob- 
tained without any need for PopIII stars. The bulk of the 
ionizing photons in those models is produced by halos with 
virial temperatures just above 10 4 K, with increasingly bet- 
ter solution if normal, PopII stars are allowed to form in 
minihalos. As shown in Fig. 1101 where we plot the UV mag- 
nitude as function of the total mass of z=7 galaxies, sources 
detectable with JWST have halo masses M h > 10 9 M Q . This 
corresponds to circular velocities v c > 50 km s . These ob- 
jects are large enough that suppression by UVB photoion- 
ization filtering is at best marginal, if not negligible at all, 
as most of the works above agree upon. Internal mechanical 
feedback might be indeed more important. This process is 
however is already included at best in the simulations when 
computing the star formation rate of individual galaxies, 
modulo the many uncertainties that still plague our under- 
standing of such phenomenon. 



9 SUMMARY 

By using high resolution simulations specifically crafted to 
include the relevant physics of galaxy formation, along with 
a novel treatment of the metal dispersion that allows us to 
follow the PopIII-PopII transition as dictated by the criti- 
cal metallicity scenario, we have been able to reproduce the 
observed UV LFs over a wide redshift range, 5 < z < 10. 
We have also shown, by combining the simulation outputs 
with a dust model previously developed for LAEs (Dayal 
et al. 2010), that dust effects at z ~ 7 should be marginal, 
although this statement depends on many details (the dust 
properties and the dust distribution scale in the ISM) that 
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Figure 10. Upper panel: Mjjy vs. total (dark+baryonic) mass 
relation for simulated galaxies at z = 7; the JWST detection limit 
is shown. Lower panel: differential fractional ionizing photon rate 
distribution as a function of total galaxy mass. The last bin on the 
left shows the contribution of larger halos too rare to be caught 
in our relatively small simulation box, computed by integrating 
the observed LF down to Muy = —25. 



are poorly constrained at this time. The general picture 
that can be drawn from our investigation is broadly con- 
sistent with the available data and therefore can be used to 
make specific predictions for the JWST. It is then useful to 
schematically summarize the main findings of the present 
work: 

• The simulated high-z galaxy UV LFs match remarkably 
well with the amplitude and slope of the observed LFs in the 
redshift range 5 < z < 10. 

• The LF shifts towards fainter luminosities with increas- 
ing redshift, mimicking a pure luminosity (or density) evo- 
lution. The faint-end slope of the LF does not vary from 
z = 5 to z = 10, keeping an almost constant slope value of 
a = -2. 

• Many galaxies at z ~ 7, especially the smallest ones, are 
virtually dust-free, and none of them shows dust extinctions 
larger than E(B — V) = 0.009. This evidence allows us to 
safely neglect the effects of dust on the UV LF. 

• The stellar population of high-z galaxies shows typical 
ages in the range 100-300 Myr at z = 5 and 40-130 Myr 
at z — 7 — 8, implying that they started to form stars as 
early as z — 9.4. These objects are enriched rapidly with 
metals and galaxies identified by HST/WFC3 show metal- 
licities ~ 1/10Zq even at z = 7 — 8. The trend of decreasing 
metallicity (and increasing spread) towards low mass halos 
indicates that small galaxies are more affected by supernova 
feedback and loose a larger fraction of the heavy elements 
they produce. 

• The relation between the star formation rate and stel- 
lar mass of simulated follows an almost linear relation with 



© 0000 RAS, MNRAS 000, 000-000 



Simulating high-redshift galaxies 13 



significant spread towards the lowest masses, implying an 
almost constant specific star formation rate. 

• The bulk of the ionizing photons is produced by objects 
populating the faint-end of the LF. These galaxies are be- 
yond the capabilities of current survey, but JWST will be 
able to resolve them up to z = 7.3. 

• Massive PopIII stars continue to form essentially at all 
redshifts and in < 20% of the galaxies. However, their con- 
tribution to the total galaxy luminosity is negligible (< 5%) 
for all objects with the marginal exception of the extremely 
faint ones; their detection will be tremendously difficult even 
for the next generation of space telescopes. 

• The typical high-z galaxies closely resemble the GRB 
host population observed at lower redshifts. This fact sug- 
gests that GRBs can be used to detect and study these ob- 
jects, providing unique information about the first stages of 
structure formation. In particular, they can be used to ex- 
tend the study of the mass-metallicity and its evolution to 
very high redshifts. Moreover, if PopIII stars end their lives 
in a GRB explosion, the detection of one of these bursts 
might represent the most promising way to directly detect 
the very first stars. 
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